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' An isotropic anti-ferromagnetic quantum state on a square lattice is characterized by symmetry arguments 

only. By construction, this quantum state is the result of an underlying valence bond structure without breaking 
any symmetry in the lattice or spin spaces. A detailed analysis of the correlations of the quantum state is 
given (using a mapping to a 2D classical statistical model and methods in field theory like mapping to the 
non-linear sigma model or bosonization techniques) as well as the results of numerical treatments (regarding 
exact diagonalization and variational methods). Finally, the physical relevance of the model is motivated. A 
comparison of the model to known anti-ferromagnetic Mott-Hubbard insulators is given by means of the two- 
point equal-time correlation function obtained i) numerically from the suggested state and ii) experimentally 
from neutron scattering on cuprates in the anti-ferromagnetic insulator phase. 
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I. INTRODUCTION 



■ The description of strongly correlated systems is one of the fundamental topics in condensed matter physics. These strongly 
^ ' correlated systems are studied mainly using low dimensional instances, where quantum phenomena are relevant and strong and 
thus one cannot expect to get an accurate understanding using classical concepts. Quantum phenomena contribute substantially 
in the Mott-Hubbard models, where the focus is on questions concerning quantum critical phenomena, metal-insulator transi- 
tions and quantum anti-ferromagnetism - with obvious connections to quantum technology. Very helpful for the description of 
any physical system is the exploitation of symmetries. Using the symmetries, simplified models (toy models) can be built to 
, understand some of the essentials of a given system. 

' In quantum anti-ferromagnetism, the traditional approach is linear spin wave theory (LSWT). Here, the assumption is that 
I the quantum system is well described by taking the classical ground state (the Neel state) and a subsequent quantization of 
T— I ■ the (linearized) deviations from that state. In the last decades, new models of anti-ferromagnets have been studied to describe 
' , particular features of Mott systems where strong interaction effects cause a strong deviation from the picture of the Neel state. 
Among these new models, there are the quantum spin liquids, which almost appear to be the opposite of the classical Neel order. 
00 : These spin liquids have at least two main characteristics that define them: (i) a spin liquid is a quantum state without magnetic 
long range order, and Neel order implies, in fact, magnetic long range order; (ii) a spin liquid is a quantum state without any 
' spontaneously broken symmetry, and the Neel state breaks translational invariance and SU{2) spin invariance. 

These models are not just an academic exercise, they try to describe states of physical systems that can appear in Nature, like 
for instance in layers of Cu02. Some of the appealing properties of these systems come from the "normal" state of the cuprate 
which is known lH 01 to be (i) a spin singlet i.e., it does not break SU (2) invariance and it does not have magnetic long range 
" ■ ' order; (ii) rotationally invariant, with a dj.2_y2 orbital symmetry; (iii) a time invariant quantum state. 

In the next sections, we will analyze a class of quantum states which have a local tensor description that enforces these 
important physical symmetries, the so called multipartite valence bond states {MVBS). One of the most frequently studied states 
with such a local tensor structure is the ground state of the Affleck- Kennedy-Lieb-Tasaki (AKLT) model lHH], a spin-1 chain 
with a Heisenberg-like Hamiltonian. The purpose of this work is to extend the AKLT construction to higher dimensions and study 
the physical properties that emerge. The extension of the AKLT model that will be studied has been described and motivated in 
detail in a prior publication, see Ref. [jBl). 

The paper is organized as follows: in Sec. 2, we describe the tensor structure underlying the construction of the ground state 
of the multipartite valence bond model. In that section, we introduce the notation and give a summary of the ideas behind the 
construction. The following sections are then dedicated to the study of the physical properties of the theoretically motivated 
model. To do so, we employ several different approaches and work out in detail what we infer and what relevance the results 
add to the theoretical model for the description of real-world physics. Along these lines, in Sec. 3, we use a mapping to a two 



dimensional statistical model to characterize the properties of the state. This mapping will allow us to extract the properties of 
the two dimensional quantum states from a one dimensional quantum system. In Sec. 4, 5 and 6, we give a detailed study of the 
correlations in the system applying well-known analytic tools (mapping to the non-linear sigma model, bosonization techniques) 
and numerical methods (exact diagonalization and variational methods). The main aim of these two sections is to try to discern 
between an algebraic or exponential decay of the equal time correlations in the multipartite valence bond state. 

Finally, in Sec. 7, we give a physical motivation of the real-world relevance of our theoretical model. We compare the two- 
point equal-time correlation function obtained from theoretical predictions of our anti-ferromagnetic model and linear spin wave 
theory with data from neutron scattering experiments of cuprates in the anti-ferromagnet insulator phase. 



II. TWO-DIMENSIONAL ANTI-FERROMAGNETIC SYSTEM. THE MODEL 

In this section, we characterize a quantum state defined on a 2D square lattice with symmetry arguments only. The aim of this 
section is to give an understandable summary of steps for the definition of the multipartite valence bond state but we will not try 
to give all the detailed ideas that brought us to the construction of the state, for that purpose we refer to iH. 

The two dimensional multipartite valence bond model is defined with the following requirements: 

• It is a real singlet of SU{2), i.e., it does not break SU{2) invariance, it does not have magnetic long range order and it is 
a time invariant quantum state. 

• It is a homogeneous, translationally and rotationally invariant state, i.e., it is invariant under any point group symmetry of 
the lattice. 

• Its local degrees of freedom are characterized by superpositions of singlet (spin-0) and triplet (spin-1) representations of 
SU{2). 

• It is the ground state of a local Heisenberg-like Hamiltonian. 

Many of the features and structures that we will find in this work can be seen as a generahzation of the AKLT model IHSl- 
The basic building block of this ground state is given by a structure known as the valence bond. We will hence call the resulting 
state the multipartite valence bond state (MPVBS). 
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Figure 11. 1 : (color online) Valence bond ground state. The first line with dark blue spots represents the physical state of the spin chain, while 
the second line represent the implementation of the state with two ancillae systems per site and a maximally entangled state between neighbor 
sites. 

Following the original works iHQ], the AKLT model describes a translationally invariant anti-ferromagnet spin-1 chain. Its 
construction is done in two steps: 

1 . In the ancilla picture, every other pair of contiguous spin-i sites is projected into a singlet state, i.e., a maximally entangled 
state with zero total spin angular momentum (valence bond). 

2. The remaining contiguous pairs are projected into the triplet subspace of two ancillary spin-i subsystems, corresponding 
to a local spin-1 site of the "physical system". 

Denoting the state of a spin-i ancilla subsystem by \a) e C^, the first step in the construction of the VBS is equivalent to fixing 
the state between neighboring ancilla spins to 

|0)= = I Ti) - UT), (ii.i) 

{a,/3}={Ta} 
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with e^i = — e^i = 1 and e|| = —e-n = 0. The projection of two spin-i subsystems into the triplet subspace is imposed by the 
definition of the physical particle to be 

(V2\ + l) a = /3=r 
(\am + ma))=l\0) a^p (11.2) 

^ [V2|-l) Q = /J=i. 

Following the valence bond construction of AKLT lUlHl, we will characterize a quantum system in a two-dimensional lattice 
lit], but in a way that is complementary to their extension of the model to two dimensions. We are looking for a uniform, 
translationally invariant singlet state in the square lattice. For this purpose, we restrict our attention to states of the following 
form (for graphical description see Fig. III. 21 and Fig. III. 3b : 



Figure II.2: The multipartite valence bond ground state. On the left we find the lattice structure describing the model, where the local physical 
degrees of freedom are located at every link. The interactions take place around every vertex and involve the four nearest neighbor sites. On 
the right, the representation of the multipartite valence bond state in the ancillae picture. 



1. On every link, the local Hilbert space is defined by a (weighted) projection of two ancillae spin-1/2 subsystems into the 
triplet and singlet representation of SU{2). 

2. On every vertex, four ancillae subsystems are projected into a singlet of SU (2). 
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Figure II. 3: Representation of the physical Hilbert space, defined at every link by the projection of two spin-1/2 ancillae subsystems, and at 
every vertex by the projection of four spin-1/2 ancillae subsystems. In the analyzed example, a singlet state is chosen at the vertex (It turns out 
to be a special case of the classical six vertex model). 



Mathematically, both conditions are imposed using a set of appropriate tensors - as in the case of the AKLT-model. First, at 
every link, the local quantum states are characterized by 



|^(A))=a(0)a°|0)+ ^ a{s)a^ \s) ^ ^L^^a'\0) + ^^L__ ^ |,), (nj) 

s={x,y,z} s={x,y,z} 

where the state |0) represents the singlet state and {|a:), \y), \z)} a basis of the triplet sector. The scalars a(0), {a(s)} give the 
amplitude of probability to find the state in the singlet and triplet sectors, respectively, such that: a{x) ~ a{y) = a{z) are the 
probability amplitudes of the triplet sector; and a" correspond to the usual Pauli matrices. 

Second, looking back at the construction of AKLT, they built a real singlet state out of real singlet valence bonds, mathemat- 
ically, they linked every second ancillary spin-1/2 subsystem with a tensor F^^^ = e^^^ that corresponds to the Levi-Civita 
tensor and is uniquely given by the symmetries of the problem. In our case, at the vertices, four spin-1/2 meet to form 
a singlet state of SU{2). Decomposing the Hilbert space spanned by these spins into invariant subspaces of SU{2), i.e.. 
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= 0+©0_©l©2, we find that the singlet sector is not unique but two-dimensional. If we rewrite the 
state at the vertex in a symmetric gauge, we find that any real singlet state can be described by 



where the angle (p determines a weight of the positive and the negative chiral singlet sectors j^. There are two special points 
(j> = {0, -j} where this tensor describes a rotationally invariant real singlet. 

It is straightforward to realize that the multipartite entangled state at every vertex corresponds to a special case of the classical 
six vertex model (see for instance 13,1^ and references therein), so that the state is decomposed in a classical structure (scaffold- 
ing) defined at the vertex and the link where we place the singlet-triplet degrees of freedom. The six vertex model characterized 
by the tensor Fj^^^^^ [(f)] lies on the critical line with anisotropic parameter A = 1. This critical line describes the phase transition 
between the disordered phase and the ferroelectric phase of the six vertex model. 



III. GROUND STATE PROPERTIES AND CORRELATIONS. EXPECTATION VALUES 



In this section, we characterize the quantum state we have just defined and obtain its correlations. Any expectation value can 
be calculated via a mapping of the 2D multipartite valence bond state to a 2D classical statistical model, and from there to a ID 
quantum mechanical model. Hence, we obtain the properties of a two dimensional quantum model from the behavior of a one 
dimensional quantum system. In the following, we give the detailed steps for this mapping. 

Given the 2D multipartite valence bond state, we map the (local) projectors to classical vertex matrices R. These matrices 
are hence derived directly from the tensors defining the projections at the Unks and at the vertices. E.g., in the calculation of 
the vacuum-to-vacuum expectation value, the projectors defining the valence bond model are mapped to two different classical 
structures: 

• At every link, a transfer matrix is defined hy E ^ J2s={o x y z} ('^'')* ® '^^ 

• At every vertex, a vertex matrix is defined hy V = T* (g) T 

and following from these definitions, the resulting i?-matrix is given by 

KZ= E ^n,n^^s,sV^:j:^n.^.^u,^. (in.l) 

The vacuum-to-vacuum expectation value in the valence bond state now corresponds to the calculation of the classical partition 
function of the derived classical vertex model; 

{m= E U^i^^^D (III.2) 

all configurations lattice 

where we sum over aU possible configurations compatible with the weights R. 
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Figure III.l: Vertex model in a classical context. The top part represents the vertex the Boltzmann weight of the vertex model. The bottom 
part represents the row-to-row transfer matrix for the vertex model. 



The classical i?-matrix or Boltzmann weight now depends on two parameters: (j), the angle in the vertex state, and A, the 
difference in the amplitude of probability of being in the triplet and singlet sector at every link. Because we would like to 
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interpret the i?-matrix as the Bohzmann weight of a hermitian Hamiltonian, we have to look at the eigenvalues of this i?-matrix 
and check for which values of the parameters they are positive definite. We restrict our attention, without loss of generality, to 
values of A > 0. The eigenvalues of the i?-matrix that can be negative are 



ri = ^ (^/3 sin (20) - 3 cos (20) 
r2 = \- (V3sin(20) - 3 cos (20) 



(III.3) 



In both cases, the i?-matrix is positive definite if G [-1, ^]. Thus we restrict the study of the vertex model to the parameter 
regimen defined by A G [0, 1] and e [f , 3^]- 

From this point, we proceed with another well-known mapping, the equivalence between D-dimensional classical statistical 
models and d = D — 1 dimensional quantum mechanical problems. In quantum mechanical problems, time evolution is carried 
out by the row-to-row transfer matrix T, with matrix elements Tab given by 



Tal 



pCiOi DC2fc2 . . . pC-ntln 

a,lC2 a.2C3 ' ' ' ^a„ci ■ 



(III.4) 



If we understand this time evolution operator as the exponential of a quantum Hamiltonian Hid = ~ logT, then the classical 
partition function Z2D is just given by the trace of the density matrix of a Gibbs ensemble Z2D = Tr ( e~^^i 
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Figure III.2; Ladder model of spin-1/2 particles. This one dimensional quantum model describes the correlations of the two-dimensional 
multipartite quantum valence bond state. For more details, please, see the main text. 



To write down the one dimensional Hamiltonian explicitly, we first realize that the indices of the i?-matrix have rank four The 
corresponding Hilbert space can be represented by two quantum spin-1/2 systems. The one dimensional quantum Hamiltonian 
can thus be described as a ladder of coupled spin-1/2 particles. To build the quantum ladder Hamiltonian we will use the spin 
operators 

and we will use a notation such that the spins in the upper row are labelled with consecutive odd integers and the spins in the 
lower row are labelled with consecutive even integers. In summary, the ladder Hamiltonian can be written as 



H^l3i ((S'l • 5^2) + (5^3 • 5^4)) + f32 [{Si ■ S3) + {S2 ■ S4)) + 133 [{Si ■ S4) + {S2 ■ S3)) + 

/34 [{Si ■ 52) (5^3 • 5^4)) + 135 [{Si ■ S3){S2 ■ 54)) + 136 [{Si ■ S4){S2 ■ S3) 



(III.6) 

2 • ^3) j , 

where the interaction involves the nearest four spin-1/2 in a plaquette of the spin ladder. 

The values of the constants f3i depend on two parameters A and 0. Initially, we will fix = ^, for two reasons: (i) for 
simplicity, in a preliminary study of the phase diagram over {A , 0} and due to the fact that these values describe a rotationally 
invariant real singlet state, (ii) as we mentioned at the beginning of this paper, the different values of describe different 
instances of the six vertex model, nonetheless all of them correspond to the anisotropy parameter A = 1 and lie on the critical 
line of the six vertex model, (iii) Our numerical study of the phase diagram {A and 0} with the DMRG method jgifioll shows 
that the parameter is irrelevant (see Fig. |III.3l l. 
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Figure III. 3: Minimum gap of the ladder Hamiltonian that encodes the behaviour of the correlation length of the 2D multipartite valence bond 
state in A and This phase diagram is built with 100 sample points using a DMRG method. Dark color describes states with small energy 
gaps and light color bigger ones. 



Then, the constants Pi{A,^) are given by 

ft 



12V9 + 6A2 + 49A4 



9 + 6A2 + 49 A4 log 9 + 7A^ log 36 - log 46656 - 2 (3 - TA^ + 4^9 + 6A2 + 49A4) log A j + 
((-3 + 7A2) log (3 + 7A2 - V9T6A2T49A4) + 3(3 - 7A^) log (3 + TA^ + ^9 + 6A2 + 49A4))] 
2(V9"+ 6A2 + 49A4 log A + (-3 + A(12 + 7A)) log ( 3+7A^+V9+6A^+ 49A1)) - 

V9 + 6A2 + 49A4 

2(V9 + 6A2 + 49A41ogA + (-3 + A(-12 + 7A)) log (S+ZAf+^ISAE+llAl)) ^ 



/32 - ^ (^log(6561/A4) - -^^^^ — ^ ' ' ■ ' gA 



/33 = T7T |log(9/A*) 

/34 



V9 + 6A2 + 49A4 

1 



3V9 + 6A2 + 49A4 




49 A4 log 9 - 7A^ log 36 + log 46656 - 2(-3 + 7A2 + 2 V9 + 6A2+49A4) log A 



- \/9 + 6A2 + 49 A4) + 3(-3 + 7A2) log (3 + 7A2 + VirM3A2 + 49A4) 

6A2 + 49A4 1og81 - 2(-3 + A(12 + 7A)) log ( 3+7A^+V9+6A-+ i9Al)) \ 
V9 + 6A2 + 49A4 / 



(111.7) 



2(V9 + 6A2 + 49Flog A + (-3 + A(-12 + 7A)) bg ( 3+7A-+V9+6A^+49A4 )) \ 

V9 + 6A2 + 49A4 J 

It is illuminating to show a plot of the values of these coupling constants as a function of the parameter A (see Fig. IIII.4l and 
Illl.St . From this first figure, Fig. IIII.4I we can deduce that the coupling between diagonal spins in the plaquette is negligible; 
and when A ^ 1 the ladder decouples into two independent Heisenberg spin-1/2 chains, since /32 — * 1 and f3i 0. This fact 
allows us to do perturbation theory around two weakly coupled Heisenberg chains [1 lj |. 

From the second picture. Fig. IIII.5I we can deduce that the values of the three coupling constants of the four spins in a 
plaquette are of the same order, i.e., none of them are negligible and there are regions where these couplings can be negative so 
that they cause ferromagnetic interaction. 
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Figure 111.4: Values of the coupling constant that involve two-body interactions. Blue corresponds to the coupling /3i along the rung of the 
ladder; Red corresponds to the coupling 02 parallel to the spin-1/2 chains; Orange corresponds to the coupling /Js between diagonal spins in 
the plaque tte. 



FOUR BODY COUPLINGS. 

COUPLINGS 




LAMBDA -A- 



Figure III. 5: Values of coupling constant that involve four body interactions. Blue corresponds to the coupling /34 along the rung of the ladder; 
Red corresponds to the coupling Ps parallel to the spin-1/2 chains; Orange corresponds to the coupling l3e between diagonal spins in the 
plaquette. 



IV. SEMICLASSICAL ANALYSIS: MAPPING TO THE NON-LINEAR SIGMA MODEL. 

We have seen that the expectation value of any operator in the 2D multipartite valence bond state can be mapped to a one 
dimensional quantum problem that is equivalent to a quantum ladder with two body and four body interactions. To g et some 
insight into the behaviour of the system, we are going to use a semiclassical analysis using the Haldane mapping lll2l [l3i [Till 
for the ladder problem ifTsIl (see also lfl6ll ). This mapping is a commonly used tool in quasi-one dimensional quantum models 
which allows to get a first description of the excitations of these quantum models. 

We should point out that although the picture that emerges with this semiclassical study follows some of the results that we 
will get with the numerical methods, the limits of validity of this approach should be analyzed in more detail. In fact, it is 
known that the presence of plaquette interactions in a ladder system gives a rich phase diagram with different critical points 
|[i7l[il[i9l|23,|2lll. Also see ll22ll23ll24ll where this semiclassical approach is used to study two-dimensional anti-ferromagnets 
at low temperature. 

Before entering in details, we fix the notation of the spin operator in the ladder 5'^,„, with ['S'°,„, 5'^„] = 
i<5/j,iyi5m.n X^c ^"'"^'^M'"' ™ &^ch a Way that it carries three indices: a G {x,y,z} labels the direction of the spin operator; 
/i S {1,2} labels the chain in the ladder and m labels the position along the chain. 

The semiclassical analysis starts with a mapping of the spin operator at every rung to a combination of a local staggered 
magnetization operator and a fluctuating one I, 

5^^™ = + (-ir+^S<^rn; 4-^,„=0 (IV.l) 

where sis the value of the representation of the spin operator, i.e., (Sfj,m^ — s{s + l). With this definition and the commutation 
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relation of the spin operators, it is straightforward to obtain 



rm jb ] - r \ ahcic . rm b ] ■ r \ ^ a 



bc,^c . 



(IV.2) 



4s2 



1, 



In the continuum limit, (ys will become the field of the 0(3) non-linear cr-model and I will become the generator of rotations, i.e. 
Hence, the set of operators that defined our Hamiltonian are mapped to 

(^^In ' 3-211 ^ 



Sin • 'S'l^n-i-i 

Sin ■ S2.n+l 
S2n ■ "Si, ,1+1 



1 \2 

-{in) -sis+1) 



1^ - (^1)" 

-^In ■ In+l H ^ S ( /„ • iPn+l — 0n ' j ~ S'^n ' '/'n+l 



2 

(-1)" . 
2 

(-1)" 



S ( ■ 0n+l - <^n ' ] ~ S Lpn ■ 0n+l 



(IV.3) 



S /„ • (^„+i + 0n ■ In+l + S 0n-0'. 



'Jn+l 



1: 

4 

■ 'fn+l + 'fin ■ In+l + S-lfin ' 0n+l, 



Assuming that if and I are slowly varying on the scale of the lattice, we calculate the Hamiltonian in a gradient expansion with 
a lattice spacing 5. Keeping term up to O (^(ff and 0{P), only since, I effectively contain a time derivative and dropping 
constant terms, we obtain 

[Sin ■ 52, 



Si. 



n+l • *J2,n+l 



(ySin ■ Si,n+lj + ^52,1 ■ 52, n+l 
Sin ■ 5*2, „+l + ^2 



(ySin ■ 8211^ ^5i,,i+i • 52,„+l 



{S2n 



S2,n 



+ 1 



(ySin ■ "Si, n+l 
(^Sin ■ S2,n+lj (^S2n ' ^i^n+l 

Finally, summing up all the terms from the Hamiltonian, 



l{x) 



3^6^ f dif 
dx 



s^S^ f dip 
dx 



(IV.4) 



-s(s + 1) {l{x) 



-(^+l)(^(x)) £ 



s3(s + l)52 (dip 

dx 



2 2 r2 

^ = /3i + /32 - s(s + l)(/34 + k)\ ^='—ip2-P3- sis + l)(/35 + f3e)) 



(IV.5) 



This effective Hamiltonian corresponds to the 0(3) non-linear cr-model with zero 0-term; it describes at every site of the 
lattice a particle moving on a sphere ((^(x))^ = 1 with angular momentum l{x) ll25l l26ll . The angular momentum takes all 

possible integer values, l{x) = 0, 1, • • • , 00. In the limit ,g ^ 1 the kinetic term ^ (ji^fj dominates over the potential term 

Jg^ (ff ) ' ^^'^ ground state is obtained choosing the smallest possible value of l{x) at every site. The first excited states 
are obtained by choosing the irreducible representation ^ = 1 at one site and / = in the rest of the chain. Since this can be done 
at every site, there is a huge degeneracy, which is broken by the potential term. It delocalizes the I — 1 excitations. The 3iV 
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degenerate first exited states become a band of / = 1 magnons, separated from the ground state by a gap, which can be computed 



using perturbation theory. In li25ll this gap was computed up to 6th order in 4y. The first three terms read 

In the weak coupling limit, a perturbative RG analysis ll27ll shows that the gap vanishes exponentially as 

A ~ -e-^'r/f (IV.7) 
9 
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Figure IV. 1: Values of effective coupling constant that define the 0(3) non-linear (T-model and minimum gap obtained from perturbation 
theory. Blue corresponds to the coupling constant g; Red corresponds to the effective velocity v; Orange corresponds to the minimum gap A 
in the strong coupling limit and green one to the gap in the weak coupling limit. 



V. BOSONIZATION OF A SPIN LADDER 



In this section, we wi ll ap ply another well-known tool from field theory, the bosonization of the spin degrees of freedom 
(see for instance ifTH [T4II28I l29tl and reference therein). This technique is complementary to the mapping to the sigma model, 
in the sense, that this method becomes exact when the coupling between spins in a rung of the ladder is weak while the one 
applied in the last section corresponds to the strong coupling limit. Nonetheless, we will see that the phase diagram covered 
by bosonization is in some sense richer than with the sigma model, giving the possibility to characterize different phases and 
to understand the transition between them. Also, it should be mentioned that bosonization techniques will allow us to have a 
qualitative description of the behavior of the system but several non-universal, i.e. model dependent, quantities will appear in 
the calculations that should be matched with numerical techniques. 

In the next lines, we will apply a method known as abelian bosonization keeping just the most relevant terms in the interactions. 
We will see that this procedure is good enough to give the correct physics that appears in spin ladders. Nonetheless, a more 
careful and detailed study of the ladder is probably needed. Several features have not been analyzed like the importance of 
marginal operators or the SU{2) invariance is not explicitly kept in our derivation; for that propose, we suppose that more 
advance methods like nonabelian bosonization would be required (see for example jSOll ). 

In what follows, we would like to look at the physics of two leg Heisenberg spin-1/2 ladders. We will start studying the 
properties of two decoupled Heisenberg chains and then we will introduce a weak coupling between the chains along the rung 
and the plaquettes of the ladder. Hence, a one-dimensional Heisenberg spin chain is defined by the Hamiltonian, 

H ~ SnSn+l = {Sn^n+l + •S'n'S'il+l) + S^S^^i ~ H^y + Hz (V.l) 

n n n 

with the usual spin operators, [S"", S^] = i6n,m£"'^'^ S^, and {a, (3, 7} € {x, y, z}. In the previous Hamiltonian, we have split 
the contribution that comes from the interaction in the xy-plane from the interaction in the 2;-spin direction, for a reason that will 
become immediately clear. 

The first step to describe a spin chain in terms of bosons is to write the spin operators in terms of Jordan- Wigner fermions, i.e. 

a, = ( n 2^rn ) - iSf) ; {ai,a^} = 5ura {a,, a™} = 0. (V.2) 

\m<l J 
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With this transformation the Hamiltonian is recast into 



From this equation we reahze that the first part of the Hamihonian, the one coming from the a::?;-interaction, is bihnear in fermion 
operators; hence, with a Fourier transformation, we can diagonahze Hxy Moreover, taking the continuum Hmit, keeping the 
lattice spacing a, and linearizing the Hamiltonian around the Fermi points (kp = the first part of the Hamiltonian has 

two independent contributions due to left and right moving fermions, 

H^y ^ ia / dx (^a^]^{x)do,aLix) - a^j^{x)d:,dR{x)j (V.4) 



with 

a{x) = ^ = e'*^^"ai^(a;) + e-''=^"at(x) (V.5) 

Now, the crucial point made by bosonization is that the low energy description of this model is given by particle-hole excitations, 
i.e., it is characterized by excitations written in terms of bilinear fermionic operators, like the electron density, which fulfill 
bosonic commutation relations 

Pr{x) =: a'i^{x)aR{x) :, Pl{x) a{{x)aL{x) :, 

iS ^ (V-6) 



where :: means normal ordering with respect to the ground state of the Hxy Hamiltonian. With these definitions (and up to 
operators that modify the expectation values with logarithmic corrections ISlfl ). the Heisenberg Hamiltonian can be recast into, 

H^H,y+H, = ^ I dx Q {pl{x) + pR{x)f + g {pl{x) - PB{x)f^ (V.7) 



with u = ^ ^(tt — l)(7r + 3), which corresponds to a renormaUzed Fermi velocity that always can be fixed to u = 1 and 
.9 = \J fq:^' this value corresponds to a weak coupling calculation, comparing the field theory methods and Bethe ansatz results, 
it is known that in the (isotropic) Heisenberg model (7=5- 

With these transformations, the fermionic operators have a simple form in a bosonic language. 



1 / 1 \ 1 / 1 \ ^^'^^ 

Pl{x) = ^ (V5 + ^ j dM^) + (V5 - ^ j dM^) 

where the c/)^ fields fulfill the following commutation relations, 

dycj),{y)] = 2m5{x - y)6^,, (V.9) 

and the Hamiltonian is recast into, 

H = H,y + H, = ^ I dx + {d,Q{x)f) (V.IO) 

with the bosonic field ^{x) = ((/>i(x) + 4'r{x)) and the dual one Q{x) = [(j)L{x) — 4>r{x)). 

Extending the relations between the bosonic and fermionic fields to the spin operators, the following equations can be found 

S'ix) = ^d^'^ix) + sin iVgHx)) 

2" ™ (V.ll) 



S+ia 



V2^ "''^ [2^ ) 
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In the next lines, we will bosonize the two legs spi n ladder, following the previous preliminary steps. Although, there is a huge 
literature on this subject (see for example ll32ll33[[34ll ), we will follow the methods used by Tsvelik and coauthors lfl7l[l8[[l9l[35ll 
. In particular, we will see that perturbed field theory for the Heisenberg ladder has only one relevant perturbation which preserve 
the SU(2) symmetry and the leading SU(2)-invariant perturbation is of dimension 1. 

For that propose, we introduce a new index a E {1,2} for every spin operator which specifiy the chain of the ladder. The 
model, that we consider, is given by the Hamiltonian 



H — -ffrung + ^leg + -Hplaquette 

— Pi (^S„ l ■ SVi,2^ + P2 (Sn^a ■ "5,1+1,0-^ + ^4 (^S„ A ' Sn 2^ {^n+l^l ' '5'„+l,2^ 

n n,a n 

n+1,1 ] iSn.2 ■ Sn+1,2] + 06 (Sn,! ' SVi+l,2j ( '5Vi,2 ' SVt+l,! 



(V.12) 



The first term in the Hamiltonian is the Heisenberg interactions between the nearest neighbor spins in a rung of the ladder, the 
second term is the Heisenberg interaction between the nearest neighbor spins in a chain of the ladder; finally, the third, fourth 
and fifth interactions describe four-body spin interactions within in a plaquette in the ladder 

The second term in the Hamiltonian is the one that we have just obtained and its bosonized expression reads 



ffieg = ^ E / {{d.<fAx)f + {dMx)y 



(V.13) 



The first term in the ladder Hamiltonian, using the definitions of the spin operators in terms of bosonic field and keeping the 
most relevant terms, can be recast into. 



J3i_ 



dx 



;(V5('i>i(a;) - $2(2:))) - cos{^{^,{x) + ci>2(.T))) + 2( 



9i(x)-92(a:) 
2^5 



(V.14) 



For the plaquette interactions, it can be seen that can be mapped to a similar bosonic interaction with an effective four body (3sS, 
that depends on the values of (34, and /Jg and non universal coefficients like the lattice spacing a, hence, it is actual value 
should be fixed comparing the numerical results and field theory methods ll20ll2lll . 



-^plaquette 



/?eff 
27rQ! 



dx 



cos i^i{x) - $2(a;))) + cos (^5 + ^2(3;))) 



Introducing the linear combinations of the fields $1 and $2^ 



$1 ± $2 



(V.15) 



(V.16) 



the decouples from the $_ contributions and the total Hamiltonian can be seen as the sum of two independent Heisenberg 
chains. 



^ I dx + {d,e+ix)f) + 



/3eff - Pi 

27ra 

P,s + Pi 



cos ($+) 



H. = ^ f dx + (9.e„ {x)f) + cos ($-) + ^ cos (e_) 

At this point, we can refermionize these models introducing a spinless Dirac fermion 

1 



(V.17) 



aR,L{x) 



\/27ra 

then, the Hamiltonian i7+ recast into a free massive fermion theory. 



exp(±i(/)+;ii,L(a;)) 



(V.18) 



H+ = ivsff I dx ( a^dxOLL - a^dxCiR ) + imt / dx ( a\aR - a\aL 



(V.19) 



with mt ~ ^^2^'" Hamiltonian also appears as a free massive fermion theory but including a pairing interaction, 

1 



cr,l{x) 



exp {±i(j)+.R^L{x)) 



i/_ = TOeff / dx ( clpxCL - ctdxCR 



+ 1 



.Pl+PcS 



27ra 



j dx {c\ 



CrCl - c{cR 



iPi /~t ~t ~ ~ A 



(V.20) 
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Finally, defining the (real) Majorana fermions, 



4 



4 



at 



the total Hamiltonian is recast into the sum of four independent free massive majorana fermions 



(m e L, R) 



E 



H„ 



(V.21) 



(V.22) 



with jTis 



3/3i+Aff 
27r 



and 



H„ 



Weff 

2 



y da; {hLdxCLL - cLRdxaR) + im J dx (aLCLR) 



All 



(V.23) 



0(1) (see 



From the above formulas, one readily obtains the line where the triplet mass vanishes, i.e. nit — > 0, when 
FiglV. l[l. This transition belongs to the universality class of the critical, exactly integrable, 5=1 spin chain (Takhtajan-Babujian 
point l36l I37I1 ) with the cental charge c = 3/2. Based on the fact that the critical line is specified in terms of three massless 
majoranas, the SU (2)2 WZNW model is the universality class describing this line. The gap, which is the mass of the majorana 
triplets, opens linearly as one deviates form criticality. Owing to the SU (2) symmetry of the model no other perturbations than 
mass terms of majoranas are allowed. 



EFFECTIVE FOUR BODY COUPLING. 




0.6 



0.8 



1.0 



LAMBDA -A- 



Figure V.l: Estimation of the effective plaquette interaction in the ladder coming from the multipartite valence bond state. In the plot, it is 
presented the sum of the absolute values of four body interactions ( | /?4 1 + | /Ss | + | /^e | ) over the interaction along the rung (3i 



VI. NUMERICAL RESULTS 



In this section, we study the ladder system resulting from the mapping defined by Eq. ( IIII.6I 1 by numerical means. The 
Hamiltonian of this system depends on two parameters (j) G [tt/G, 2tt/3\ and A G [0, 1] which fix the six parameters of Eqs. 
(HILTT i. 

In the following, we present results concerning the exact diagonalization of ladders up to 28 spins or length L = lA for 
different values of A with open boundary conditions and using DMRG calculations. Several points seem to be remarkable. First, 
there is a region of criticality around the point A ^ 1; as the system approaches this region, it can be described by two decoupled 
Heisenberg chains, a critical system. Second, the first excited state is three-fold degenerate. Third, in Fig. IIII.31 we plot the 
minimum gap as a function of the parameters (0, A). From the plot, it seems that the angle </> is an irrelevant parameter at least 
for a substantial part of the phase diagram {<f>, A). Nevertheless, there are special regions that require a careful analysis. 

To achieve some knowledge about the behavior of the system in the thermodynamic al limit, we used the DMRG method 

[l3l to find the system ground state and the first excited state energy. We sample the two-dimensional parameter space with 
one hundred points, and for ladders of up to L = 60 rungs, retaining m = 40 states to compute the gap of the system (we 
checked convergence of the results till m ~ 100 for some typical points). Typical results are shown in Fig. IVI.ll These results 
have been compared with exact diagonalization calculations for small instances of the problem, e.g. L = 14. The complete 
phase diagram for L = 60, m = 40 is displayed in Fig. IIII.3I 

From the phase diagram it is clear that the gap is closing for A ^ 1. Due to the closing gap, the correlations in the ground 
state become stronger and thus the numerical simulations require more and more resources. We moved to a ladder with L = 100, 
retaining m = 60. In Fig. I VI. li the gap values for A = 1 and cj) ~ 1.15 are shown. In the inset, it can be also seen the results in 
the regime A = | and (p — ^ which signal a different behavior 
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FIRST GAP (A, 0) 

ENERGY 




0.02 0.04 0.06 0.08 LENGTH 

Figure VI. 1 : Results of DMRG calculations of first energy gap as a function of the length and the number of states kept in the simulations, with 
(A, 0) = (1, 1.15) in the main plot and (|, ^) in the inset. Blue points correspond to a number of states m = 40, purple points to m = 60 
and golden points to m = 100. In the main plot the golden points follow a linear behaviour expected in a critical theory. 



Exact diagonalisation is used here to provide exact information about the low energy spectrum of the system with the size up 
to N=28 spins and to benchmark the DMRG results. Exact diagonalization by construction is an unbiased and reliable method 
but its main drawback lies on the limitation of the systems that can be analyzed due to the exponential growth of the Hilbert 
space with the system size. The developed code utilizes the PETSc |38] and SLEPc ll39ll libraries and is run on a Blue Gene/P 
supercomputer utilising up to 2048 processing cores. 

The question that we try to solve with exact diagonalization is the behavior of the minimum gap of the system with the length 
and the perturbation A. We already got some a priori knowledge about the possible scenarios that can appear due to the field 
theoretical analysis. In the case that the system is gapless, in the thermodynamic al limit, the minimum gap should present a 
linear dependence with the system size as it is shown in the Fig. IVI.2I 



FIRST GAP(A=0.3) 




0.08 0.10 0.12 0.14 0.16 0.18 0.20 LENGTH 




0.08 0.10 0.12 0.14 0.16 0.18 0.20 LENGTH 



FIRST GAP (A=0.5) 




0.08 0.10 0.12 0.14 0.16 0.18 020 LENGTH 



FIRST GAP (A=0.9) 
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Figure VI.2: Results of exact diagonalization of first energy gap as a function of the length and A for length up to L = 14. All the plots show 
a clear linear dependence of the gap with the inverse of the length of the ladder. The lines correpond to the best linear fits and are just to guide 
the eye 



Another convincing argument about the gapless system is to show how the gap scales with the perturbation A. Now, we know 
that perturbed conformal fields theory shows that there is only one relevant perturbation which preserves the SU (2) symmetry of 
the system with scaling dimension one, which means that close to the A = 1, where we have two decoupled Heisenberg chains, 
the gap of the system should open linearly if the system remains gapped in the thermodynamical limit. In Fig. IVI.3I we have 
plotted the first energy gap as a function of A~^, there, we can see that the plots does not show a Unear dependence of the gap 
with the perturbation A. 
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MINIMUM GAP (L) 



MINIMUM GAP (L) 




Figure VI.3: Results of exact diagonalization of first energy gap as a function of the length and A. On the left side, it is shown the plots for 
system sizes L — {12, 14}. One the right side, it is plotted a closer look of the results for L = 14. The blue line correponds to the best linear 
fit, the red line corresponds to the best quadratic fit; they are just to guide the eye but it seems that the plots does not show a linear dependence 
of the gap with the perturbation A. 



VII. NEUTRON SCATTERING USING THE 2D MULTIPARTITE VALENCE BOND STATE 



This section is dedicated to the study of the real-world relevance of the theoretical model. To this aim, we compare the two- 
point equal-time correlation function obtained from the theoretical predictions of our anti-ferromagnetic model and the linear 
spin wave theory with data from neutron scattering experiments of cuprates in the anti-ferromagnet insulator phase. 

In the context of high-Tc superconductivity, the investigation of the different competing phases of a potentially superconduct- 
ing material and the transition between them is of strong interest. In the limit of low temperature and no doping, a very simple 
structure is dominant: the Neel order, the signature of an anti-ferromagnetic long-range ordered state. However, this phase is 
not the interesting one. Changing the environmental conditions or doping might cause a phase change, in some cases leading 
to superconductivity. These normal phases do not show long range order and have also other interesting structural properties. 
Among these is the formation of stripes, whose occurrence is conjectured by experimentalists to be linked to high-Tc supercon- 
ductivity ll40l l4lll42i I43II44I1 . In this section we want to investigate the multipartite valence bond state from this point of view. 
Since it is lacking long range order, it might be a candidate for the description of normal states of a wide range of materials. 
To give an example, we might be able to model cuprates, which are assembled from spin-1/2 sites, but might exhibit - in a 
low-energy sector of the normal phase - the behaviour of a substrate made of spin-1 sites. Another example is the nickelates 
which are actually assembled directly from spin-1 sites. In this case, we might be able to grasp even the full theory, not only an 
effective one. 

The investigation carried out in this section is the numerical analysis of a hypothetical neutron scattering experiment, per- 
formed on the MPVBS. It is then compared to experimental data as well as predictions of other theoretical models. Our focus 
will be on the short-range entanglement features exhibited by the experiment (and the predictions), because these features enter 
the construction of the MPVBS on a very fundamental level. 

The aim of and the essential idea behind a scattering experiment is to find out which degrees of freedom (to absorb energy 
and momentum, etc., from the incident scattered beam) a target system has, and from this to draw conclusions about the internal 
structure of the target. In principle, a theoretical model derived from first principles should provide us with these degrees of 
freedom and therefore with the possible outcomes of the interaction of beam and target. 

The central quantity for our following analysis is the experimentally and theoretically accessible differential cross section 

J ground state^any 

We assume a setting where the incoming and outgoing beams are waves of definite momentum. We are not interested in the final 
state of the target system and choose the ground state of the anti-ferromagnet as our target. The fact that we are not interested in 
the final state liberates us from the task to determine the exact Hamiltonian; when traced out, it is not important. In rotationally 
invariant states we obtain the formula 

"5 ' / ground state^any ^i.rj 

which implies a relation between the two-point correlation function and the differential cross section. For further references and 
details, we refer to e.g. ll45ll . The correlation functions are accessible in calculations involving the MPVBS. In this section and 
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using this relation, we will now make further use of these (numerically accessible) quantities and hence be able to compare them 
to cross sections of other models, techniques and to experimental data. One of these alternatives techniques is the linear spin 
wave theory. 



A. The linear spin wave approach in an anti-ferromagnet and alternative theories 

Assuming that we have akeady found the right model for the description of a solid state system, a big challenge remains: the 
problem of performing the quantum mechanical calculation in the thermodynamic limit, i.e., for objects of macroscopic size that 
are observed in the laboratory and usually consist of ^ 10^^ particles. In the context of scattering experiments, but not only 
there, a good way to reduce the complexity of the macroscopic system is to make certain a priori assumptions about its degrees 
of freedom. One of the conventional assumptions that can be made and that will be used in this paper as a reference is the idea 
that the excitations of spin waves are the essential events that take energy and momentum from the incident beam of scattered 
particles. 

If we assume that the system is properly described by the Heisenberg Hamiltonian on a square lattice 

= J ^ Sj • Sk, J > 0, 

where (j, k) is a sum over nearest neighbors, then in the corresponding classical ground state every two neighbors would be 
aligned in parallel, but with opposite orientation. This leaves us with two sub-lattices (1 and 11) with equally oriented spins, one 
sub-lattice pointing up and one down. This ordering is called the Neel ordering and the corresponding state with this ordering 
is called the Neel state. The deviations from this order can be described in terms of modes with definite wavelength. The 
quantization of these modes leads to another quantum theory of the state. The linearized version of this theory is the traditional 
description of the quantum mechanical anti-ferromagnet: linear spin wave theory. 

Hulthen ll46ll carried through the quantization of these spin waves, which subsequently were used by Anderson ll47ll in his 
approximate theory of the quantum mechanical anti-ferromagnetic ground state. They finally obtain the proportionality 

exp hik (j - 1)] cx V(l-7k)/(l+7k), 

j,l, boih from /,// 

where 7k = J2i '^'^^ i-^-' V ^ 7k) / (1 + 7k) is proportional to the differential scattering cross section into the direction of 
k. 

Another approach worth mentioning is the superposition of the classical ground state of the anti-ferromagnet - the Neel state 
- with fluctuating valence bonds. This approach is a construction of (an approximation of) a ground state that offers a qualitative 
understanding but fails to make the right predictions in certain regions of Q-space. We will not elaborate on this approach here, 
but account for its predictions later, if appropriate. 

As it has already been pointed out in ll48ll . there are many more models and techniques to make predictions beyond linear 
spin wave theory and the multipartite valence bond state. A very important technique is, of course, quantum Monte-Carlo, 
which gives the right answers (or at least being_^very close), as indicated by simulations at certain points of interest in Q- 
space which were performed and discussed in 114811 . The approximate correctness of quantum Monte-Carlo simulations for this 
model gives confidence that the basic assumptions about the model are justified and we can concentrate on finding its solutions. 
However, the quantum Monte-Carlo simulations are not suited to understand the features of the solution qualitatively and are 
hence unsatisfactory. We also stress again that the Heisenberg Hamiltonian is not the optimal choice either, since it is known that 
plaquette interactions take place in the real-world anti-ferromagnetic substrates ll49ll . However, the family of Heisenberg-like 
Hamiltonians is large and includes alternative anti-ferromagnetic Hamiltonians; also the parent Hamiltonian of the multipartite 
valence bond state belongs to this class. 



B. The experiment 

Recently, the various scattering cross sections of a two-dimensional square-lattice anti-ferromagnet have been measured ex- 
perimentally in great detail ll48ll . offering a great opportunity to test a wide range of theoretical models. 

The substance to be measured was a single crystal of CFTD [Cu(DCOO)2-4D20], with lattice constants a = 8.113, b ~ 8.119 
and c = 12.45. In an effective model of CFTD, the essential sites (i.e., those that are responsible for the properties of the 
substance) are given by the electrons in the d-shells of the copper ions, naturally being spin- 1/2 particles. In the ionic picture, the 
electrons are strongly localized and there is almost no coupling between two layers in the c-direction. Moreover, the effective 
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interactions between these sites cause the crystal to be anti-ferromagnetic. Hence CFTD is an almost perfect realization of a 
two-dimensional square lattice Heisenberg anti-ferromagnet. 

In the experiment, several assumptions were tested, (i) One assumption is the idea that the quantum anti-ferromagnetic 
ground state is essentially the classical counterpart with minor quantum corrections, which - once found - are moreover well 
understandable conceptionally. This assumption expresses the hope that the quantum anti-ferromagnetic ground state can be 
understood qualitatively, (ii) The second assumption that was tested is the idea that spin waves are the essential degrees of 
freedom of excitation in the anti-ferromagnetic ground state. 

We furthermore notice that both linear spin wave theory and the experiment assume the interaction to be described by the 
Heisenberg Hamiltonian. This assumption makes these models different from the MPVBS, because its parent Hamiltonian 
contains plaque tte terms. 



C. Comparison of theoretical predictions with experimental data 

We compare the experimental data with the predictions of the linear spin wave theory and the two-dimensional multipartite 
valence bond state. In both cases, the differential cross section is proportional to the Fourier transform of the two point function. 
The calculation of this two point function is straightforward but - without further simplifications - naturally becomes more and 
more tedious as the system size grows. In the case of the multipartite valence bond state, we did not use any further assumptions 
about the degrees of freedom but calculated the two-point function directly. In fact, we fix the parameters (A, (p) that defined 
the MPBVS to the values {\,^)- We have performed the calculations up to a grid size of 8 x 8, and we use the results of this 
lattice for the analysis. Finite size effects are an issue in such a small lattice and will be most notable for wave vectors of large 
magnitude, i.e., in the center of the Q-space, which is given by the interval (0, 2tt) x (0, 27r). 

We will now discuss the cross sections in greater detail. Initially, the cross section is given as a function 

(0,2^)x(0,2^)9Qk^ (Q). 

Va"/GS^any 

For convenience, we will follow this function on a path - being depicted in Fig. lVll.fl - that crosses several interesting regions 
of Q-space; the reasons for these regions being interesting will be explained later. This way we arrive at a two-dimensional plot 
of scattering cross section versus length on the path, which is more accessible than a three-dimensional version. The shape of 
the path is motivated by experimental work and the data provided therein. 
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Figure VII. 1 ; The path to follow. 

Our discussion starts with the raw qualitative features of the predicted and experimental cross sections using Fig. IVII.2I In this 
plot the cross sections of the linear spin wave model, the multipartite valence bond state and the experiment along the path are 
given. The units are arbitrary and all curves and data points are normalized to give the best possible match with the experiment 
overall. 

We see that LSWT predicts the experimental data with excellent quality on a crude scale. The singularity at (tt, tt) is predicted 
properly. Moreover, there are no huge qualitative differences to be seen, although the cross section is false in some regions, 
showing a slightly wrong tendency there. Most interestingly, in the first part of patch 2 of the path the cross section is generally 
too high, this feature being reflected again in the last part of patch 5, which covers a similar region. A comparison of LSWT 
with the experimental data in this region shows that close to the point (27r, 0) LSWT predicts a too small scattering cross section; 
and close to (tt, 0) the prediction is too large. The error of the prediction shows an almost linear relation along patch 5. This 
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Figure Vn.2: Cross sections along the path. 
Features of the predictions of linear spin wave model, the multipartite valence bond state and the experimental data. The inset shows the same 
data, normalized such that the value of the LSWT equals one everywhere. G. Aeppli's group provided the experimental data published in ll48h . 



feature is closely linked to the prediction of linear spin wave theory that at the point (27r, 0) the cross section should be zero. 
Considering the data, this is most likely a qualitatively wrong statement about the cross section, based on wrong assumptions 
about the degrees of freedom available and hence about the electron structure at these length scales. Generally, we observe that 
the cross section in the central region of the wave vectors, i.e., close to the point (tt, tt) are more precise in their predictions 
than those at the zone boundary, the latter belonging to wave vectors resolving features of the electron structure of small scale; 
a typical length being around the size of a small multiple of the distance Cu-Cu in the lattice. 

Like the LSWT, the MPVBS follows the experiment very well, considering the finite size effects of the simulation. Please note 
that there are two regions of the plot that we should consider separately in the context of our simulation of the MPVBS. Finite 
size effects divide the zone into two regions of different precision: like in the case of LSWT we have to distinguish between 
wave vectors in central and peripheral regions, the latter being more precise than the former numerically, since finite size effects 
have a great impact in the central region. 

The first region is represented by the patches 1, 2, 3 of the path. The singularity at (tt, tt) is not predicted. We believe that the 
deviations in this region are mostly due to the finite size of the lattice that was used to perform our calculations, and not due to 
the qualitative features of the MPVBS. However, we were unable to prove this, since direct calculations with larger lattices have 
so far not been feasible, and in principle one should go to the thermodynamic limit. The patches 4 and 5 of the path belong to 
wave vectors at or close to the boundary of the zone and hence are not strongly affected by the small lattice size. Please note 
that this is the region of the most interesting deviations of LSWT from the experiment. Here, one qualitative difference between 
LSWT and the MPVBS is apparent; at the point {271, 0) the LSWT scattering cross section vanishes completely and becomes 
comparatively small already in a comparatively large region near this point. The MPVBS shows a completely different behavior; 
the scattering cross section stays finite. Unfortunately the experiment does not provide data in this region. 

As has been pointed out before, the patches 4 and 5 of the path are most expressive, but also the region where patch 1 and 
2 meet, the finite lattice size does not have a strong impact. In the latter region, we see that LSWT predicts the experimental 
data with an error of 100%. The MPVBS, however, shows the right tendency. The predicted value is less than that of LSWT. 
The comparison of the qualitative behavior of the LSWT, the MPVBS and the experimental data is possible in an optimal way 
in a plot that shows the data from Fig. I VII. 21 normalized to the predictions of LSWT. Here the tendency in the patches 4 and 5 
is most apparent. Please notice again that the patches 1 , 2 and 3 are severely affected by the finite size of the lattice used for 
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computations. 

What is the significance of patch 5 and the meeting points of patches 1 and 2? The cross section in this region of the zone 
reflects the electron structure at short distances of the lattice. As has been conjectured already in ll48ll . the entanglement features 
on this length scale in the Heisenberg anti-ferromagnet and hence also in CFTD might be completely different from those of 
the Neel state and the spin wave excitations. On the one hand, the Neel state is a classical state with no entanglement and also 
the excitable degrees of freedom, the magnons, show the entanglement structure of collective modes of excitations of harmonic 
oscillators. 

On the other hand, the MPVBS is constructed to have the entanglement structure of (two-electron) singlets, established 
between nearest neighbors. This assumption about the electron structure obviously allows for better predictions at the small 
length scale. An assumption similar to the construction of the MPVBS is made when we start with the Neel state as the ground 
state but allow for a superposition with resonating valence bonds, in the ground state as well in its modes of excitation. As 
shown in ll48ll . this construction gives a very good match at point (tt, 0), but results at the same time in a prediction that is far 
too low - though qualitatively right - when approaching the point (27r, 0). This wave vector pointing in the diagonal direction 
of the lattice, where no valence bonds are assumed in the construction. Obviously some degrees of freedom are missing in this 
description, which again are included in the MPVBS. Please note in this context that the effective interactions in CFTD arise 
under mediation by formate groups in between the copper ions, a situation which is considered and reflected in the construction 
of the multipartite valence bond state. 



VIII. CONCLUSIONS 

In this paper we carried out several approaches for the analysis of the multi partite valence bond state described in Ref. 
This state is based on a local tensor description that enforces important physical symmetries. It can hence be interpreted as being 
a two-dimensional generalization of one of the most frequently studied states with such a local tensor structure, the ground state 
of the Affleck-Kennedy-Li eb-Tasaki (AKLT) model iHQl, a spin-1 chain with a Heisenberg-like Hamiltonian. Moreover, it can 
be seen as a complementary generalization to their extension to higher dimensions. 

Our main goal in this paper was to analyze the behavior of the two point correlation functions of the multipartite valence bond 
state model and to compare it with the classical linear spin wave theory (LSW) and experimental data. 

To achieve this goal, we proceeded in two directions: first, we applied a mapping to a quantum ladder problem, in order to 
discern between an algebraic and an exponential decay of the correlations. In this step, we use analytical tools (mapping to 
the non-linear sigma model and bosonization techniques) and numerical tools (exact diagonalization and variational methods). 
Another strategy that should be investigated in the future, and we leave open, it is the study of the equivalent two dimensional 
classical vertex model, its integrability properties and the relations with the Yang-Baxter equations. 

Second, we calculated the equal time two point functions and compared it with both the results of LSW and data from a 
neutron scattering experiment on a Mott-Hubbard anti-ferromagnet insulator material. 

In the context of the description of two-dimensional solid states in the anti-ferromagnet insulator phase, our model can be 
seen as an alternative to the linear spin wave approach. In fact, the scattering cross sections of both models differ qualitatively 
in important regions. These regions correspond to structures that the LSW approximation fails to describe properly and that are 
conjectured to correspond to important qualitative features of the electron structure in the actual solid. More precisely, these 
features are believed to be related to short-range entanglement structures of the electrons in the solid and hence reflect a design 
goal of the generalized AKLT model that we give. 

Moreover, the multi-partite valence bond state offers a new and unusual constructive approach for toy models in solid state 
physics, providing unusual points of view on the structure of high-Tc superconductivity. 
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